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Organic semiconductor density of states controls 
the energy level alignment at electrode interfaces 

Martin Oehzelt^'^, Norbert Koch^'^ & Georg Heimel^ 



Minimizing charge carrier injection barriers and extraction losses at interfaces between 
organic semiconductors and metallic electrodes is critical for optimizing the performance of 
organic (opto-) electronic devices. Here, we implement a detailed electrostatic model, cap- 
able of reproducing the alignment between the electrode Fermi energy and the transport 
states in the organic semiconductor both qualitatively and quantitatively. Covering the 
full phenomenological range of interfacial energy level alignment regimes within a single, 
consistent framework and continuously connecting the limiting cases described by previously 
proposed models allows us to resolve conflicting views in the literature. Our results highlight 
the density of states in the organic semiconductor as a key factor. Its shape and, in particular, 
the energy distribution of electronic states tailing into the fundamental gap is found to 
determine both the minimum value of practically achievable injection barriers as well as their 
spatial profile, ranging from abrupt interface dipoles to extended band-bending regions. 
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Organic electronics has tremendous potential for extending 
and complementing the range of conventional semicon- 
ductor applications and has, to some extent, already 
entered the market. This is especially true for organic light- 
emitting diodes (OLEDs)^"^ in displays of mobile devices 
and, very recently, in television panels. In the fields of solid 
state lighting, integrated circuits based on organic field-effect 
transistors (OFETs)^"^^ and organic photovoltaic cells 
(OPVCs)^^"^^, commercial breakthrough is still ahead, with 
considerable effort in fundamental as well as applied research 
currently undertaken. All these devices have in common one or 
more organic semiconductor layer (s), possessing (as every 
semiconductor) a non-vanishing energy gap between occupied 
and empty electronic states. Electrical connection of the organic 
semiconductor(s) to peripheral circuitry, however, must 
inevitably be realized by metallic contacts, with a wide range of 
elemental metals and degenerately doped inorganic or organic 
semiconductors employed today^'^^'^^. These are characterized by 
a continuum of electronic states filled up to a common electron 
chemical potential, often referred to as the Fermi level, E^. 
Depending on the intended function of the organic electronic 
device, charges must either be injected from the metallic electrode 
into the organic semiconductor (OLED/OFET) or they must be 
collected from the organic semiconductor by the metallic 
electrode (OPVC)^^'^^. 

In Fig. la, a typical situation for the energetic alignment of all 
relevant electronic states is sketched schematically. We immedi- 
ately realize that, in such a situation, charge carriers encounter 
energy barriers, Ae for electrons and A^ for holes, if they are to be 
injected and that, for their extraction, corresponding energy 
losses occur. This has led the community to acknowledge that 
minimizing these barriers/losses and, in a broader sense, 
controlling interfacial energy level alignment in general, is a key 
factor for improving device performance 

If, for a specific material pair, the situation sketched in Fig. la 
is found unsatisfactory, the obvious strategy is to reduce barriers/ 



losses by either moving towards the occupied or unoccupied 
molecular levels (for example, by modif)^ing the electrode work 
function or by choosing a different electrode material) or to move 
one of the molecular levels closer to E^ (typically by choosing a 
different organic semiconductor). 

As the ultimate limit of such a strategy, one could assume that 
it was possible to eliminate barriers/losses altogether (Fig. lb). In 
practice, however, this was only observed for monolayers of select 
molecules^^"^"*, but could never be achieved for multilayer 
organic semiconductor films of device -relevant thickness. There, 
in contrast, non-vanishing energy barriers were found to remain 
(Fig. Ic). Several models have been developed in the past to 
explain this phenomenon^^"^^'^^"'*^ each of which will be 
discussed to some detail in corresponding sections below. They 
differ not only in their intended range of applicability (from a 
single organic semiconductor layer^^ to thick films^^'^ '^^'^^) but 
also in the envisioned regimes of electronic coupling strength 
between organic and electrode (from weak coupling with a 
passivated electrode^^'^^'^^'^^ to intermediate coupling with an 
atomically clean metal^^) as well as in the assumed energy 
distribution of the electronic states in the organic (from 
discrete^^'^^'^^ to continuous^^'^^). 

Here, we present results from a numerical electrostatic model, 
which allows treating all these various scenarios on equal footing. 
The full wealth of all encountered interfacial phenomena is 
described qualitatively and, provided required input data are 
available, even quantitatively^"*. With previous models naturally 
emerging as different limiting cases, we consolidate diverging 
views in literature regarding the physical origin of non-vanishing 
injection barriers (Fermi level pinning by metal-induced 
density of gap states^^ versus finite-temperature occupation of 
discrete molecular levels^^'^^) as well as their spatial profile 
(interface dipole^^'^^ versus band bending^^'^^), and we provide a 
single, coherent framework for understanding the energy level 
alignment at interfaces between metallic electrodes and organic 
semiconductors. 




Metal Organic Metal Organic Metal Organic 



Figure 1 | Sketch of the energy level alignment at organic/electrode interfaces, (a) Schematic showing the Fermi-Dirac occupation function with 
Fermi level for the metallic side of the interface. The difference in binding energy (BE) between and the onsets of occupied and unoccupied density of 
states (DOS) on the organic semiconductor side are the charge injection barriers, Ah for holes and Ae for electrons, respectively, (b) Example of a possible 
initial situation for the iterative electrostatic model described in the text. The energy differences between the onsets of the DOS derived from the 
highest occupied molecular orbitals (HOMOs) and the lowest unoccupied molecular orbitals (LUMOs) to the vacuum level (VL) define the ionization 
energy (IE) and the electron affinity (EA) of the organic semiconductor, respectively. The energy difference between and VL is the work function of the 
organic film (Wffjim). (c) Final energy level alignment after convergence of the iterative scheme discussed in the text. With increasing distance z 
from the electrode, the DOS increasingly shifts to higher BE by the local electron potential energy —eV(z), thus increasing Ah and decreasing Ae with 
respect to the initial situation. 
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Results 

Introducing the model. Our model is inspired by the earlier 
work of Vasquez et al?^'^^ and, more recently, of Blakesley 
et alP'^^y as well as of Greiner et alP and Ley et al?^ Going 
beyond refs 23,25, however, we take into account a continuous 
density of states (DOS) both for the occupied and the unoccupied 
levels in the organic semiconductor. For simplicity, we start with 
the common choice^^'"*^'"*^ of one Gaussian peak, G{E — £h>^^h)> 
arising from the energy distribution of the highest occupied 
molecular orbitals (HOMOs) and one, G(E — Ei^.Gi), from the 
distribution of the lowest unoccupied molecular orbitals 
(LUMOs) as indicated in Fig. lb. To unambiguously define 
these peaks, both the energies of their maxima, and £l> and 
their standard deviations, cth and ctl, need to be specified. The 
latter can be approximated, for example, by the width of the 
HOMO-derived feature observed in ultraviolet photoelectron 
spectroscopy (UPS). Here, we adapt the constant value of 
(Ju — <^L — ^-^5 eV, found for the archetypical organic 
semiconductor pentacene"*"*. From UPS and inverse 
photoelectron spectroscopy also the ionization energy 
(IE = 4.8eV)^^ and the electron affinity (EA = 2.6eV)^5 of 
pentacene are known. These quantities are typically determined 
from the energy difference between the vacuum level (VL) and 
the low-binding energy onset of the UPS spectrum and the high- 
binding energy onset of the inverse photoelectron spectroscopy 
spectrum, respectively. Choosing the onset of a Gaussian peak to 
be the energy, where the tangent through the inflection point 
crosses the base line, puts the peak centres 2(7 away from the so 
obtained onsets, that is, at £H = 5.3eV and £L = 2.1eV, 
respectively. Also indicated in Fig. lb is the work function, 
Wffiim, given by the energy difference between the Fermi level and 
the VL outside the organic semiconductor film. Note that, in the 
specific situation shown, Wf^i^ is actually higher than the IE of 
the organic semiconductor, where electrons reside at energies 
above E^ before charge equilibration. It is, therefore, clear that in 
the course of establishing a common electron chemical potential, 
charge will flow across the interface. To quantif)^ the amount of 
transferred charge in the general case, we discretize the organic 
semiconductor film into layers of thickness Az and distance z 
from the electrode surface. We then calculate the charge density 
p{z) in each of these layers by integrating 

{+00 
J dE-fu{E)-Du[E + eV(z)] 

(1) 

+ 00 ^ ^ ^ 

- j dE-f^{E)-Di,[E + eViz)][ 

— 00 J 

with e the (positive) elementary charge, n the number of 
molecules per unit area and per discretization interval Az. As 
the IE and EA values reported above pertain to pentacene 
molecules 'standing' on the electrode with their long molecular 
axes almost parallel to the surface normal"*^'^^, we accordingly 
choose Az=1.42nm, the distance (iooi between the 
crystallographic (001) planes in the pentacene single crystal^^; 
Az then corresponds to the height of one molecular layer. For the 
sake of simplicity, we consider layer-by-layer growth of pentacene 
only, with n = 4.1 x 10^^ molecules/(m^ • Az). Note, however, that 
fractional coverages and different growth modes could easily be 
incorporated by accordingly choosing layer- specific values for n. 

The first integral in the curly brackets of equation (1) yields the 
number of holes in the occupied DOS of the organic 
semiconductor and the second integral the number of electrons 
in the unoccupied DOS. While of noticeable consequence only for 
highly degenerate systems, we, nevertheless, follow refs 23,25 by 



allowing, in contrast to refs 27,28, only single charging of each 
molecule, which slightly modifies the Fermi function for each 
individual energy level under consideration to 

1 (2) 

Here, and gi^ are the HOMO and LUMO (spin-) 
degeneracies, both 2 in the case of pentacene, and j] ^l/k^T with 
/cb the Boltzmann constant and T the temperature, set to 300 K 
throughout the paper. Finally, the energy distributions of 
HOMO- and LUMO-derived levels are Du[E + eV{z)] = 
G[E-Eu + eV{z\ Gu] and DdE + eV{z)] = G[E- E^^^ eV{z\ 
ctl] in our model case, unless otherwise noted. We stress, however, 
that (i) any peak shape can be used instead of a Gaussian, (ii) each 
discretization interval could have a different DOS to account, for 
example, for image- charge screening by the metal, (iii) more 
molecular states can easily be taken into account by adding the 
respective terms in the curly brackets and (iv) potential 
temperature effects can be included in the standard deviation 
a = (7(T); for the values used here, 2a = 0.50 eV is about 20 times 
bigger than k^T at T=300K and, therefore, the temperature 
dependence of /h/l(£) has only little impact on level alignment. 

Finally, the electrostatic potential across the organic semi- 
conductor film, V{z)y is first set to zero for calculating an initial 
guess of the charge density via equation (1). Once known, this 
p{z) serves to find the next V{z), which, going beyond refs 23,25, 
we obtain by numerically solving the generalized one- 
dimensional Poisson equation 

V[£(z)Vy(z)] = - ^ (3) 

with the discrete-difference method on a variable grid, using the 
second- order, central- difference scheme with mixed Dirichlet 
[y(0) = 0] and von Neumann [V'(d) — 0] boundary conditions, 
where d is the total thickness of the considered organic film and 
£0 the vacuum permittivity. For our first illustrative case of a 
homogenous pentacene film, we assume a static dielectric 
constant £ of ~3.6, as recently determined by impedance 
measurements'*^; note, however, that each discretization interval 
could, in principle, also be given a different dielectric constant. 
The DOS of each discretization interval is then shifted in energy 
by the respective — eV{z)y as shown in Fig. Ic, and the next p{z) 
is again generated by occupying all states in the organic up to Ep. 
This self-consistent process is iterated until the difference in 
charge density between cycles is less than one electron per layer 
and square meter and the final injection barriers are obtained as 
Ae = £F-EA and Ah^IE-Ep. 

Some illustrative results. To illustrate the results delivered by the 
method just described, we first consider 20 layers ( ~ 30 nm) of 
standing pentacene with set to lie right at the energy of the 
HOMO peak maximum. The self-consistent solution of the 
Poisson equation for this situation is depicted in Fig. 2a,b, 
showing the charge density and the electron potential energy, 
respectively. The charge density rapidly vanishes with increasing 
distance from the electrode surface and reaches almost zero at 
5nm. The corresponding potential energy, which locally shifts 
HOMO, LUMO and VL by the same amount (IE and EA remain 
constant), is characteristic for band bending— a term somewhat 
liberally borrowed from inorganic semiconductor physics despite 
the fundamental difference that, instead of actual bands there, a 
distribution of localized states is assumed here. Although this 
phenomenon has often been attributed to (impurity or 
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intentional) doping^^'"*^'^^ in the past, it has recently been argued 
to arise also in intrinsic organic semiconductors, both of the 
conjugated-polymer^^ and the small-molecule varie 
trying to relate our results to experimental observations, it is 
important to recall that the work function measured by, for 
example, UPS or Kelvin probe only reveals the electron potential 
energy —eV{z>d) just outside the organic film (highlighted in 
Fig. 2b). As pointed out in refs 50-53, however, the sequence of 
— eV{d) values obtained by increasing the total thickness d of an 
organic film in a step -wise manner and probing the surface 
potential at each step does not, in fact, recover the electron 
potential energy profile — eV{z) within the organic film for any 
thickness (Fig. 2c). 



Fermi level pinning. With this in mind, we now plot in Fig. 3 a 
the total potential energy change, — eV{d)y across pentacene films 
of increasing d^ for a series of different Fermi level energies. 
Clearly, when £p lies around mid-gap of the organic semi- 
conductor, only very little charge is transferred from/to the metal 
electrode and, consequently, the potential energy profile remains 
close to flat. Upon Ep approaching either the HOMO or the 
LUMO levels, however, pronounced charge transfer and, 
accordingly, ever higher and steeper potential energy changes 
across the organic film are obtained. Practically, the Fermi level is 



shifted by choosing electrodes of different work function, Wfeiec> 
defined by the energy difference between £p and the VL just 
above the bare electrode surface with no organic material present. 
The energy barriers/losses for charge carriers, Ag and A^, are then 
related to Wfgiec and — eV{d) through 



Ae = [Wfelec + ID - eV{d)] - EA, 

Ah = IE-[Wfelec+ID-ey(^/)], 



(4) 
(5) 

where ID is the interface dipole. Either because of contamination 
through atmospheric and/or solvent exposure or owing to 
deliberate passivation with inert buffer layers, the ID is often 
found to be (close to) zero for practically employed electrodes 
in organic electronics. However, microscopic charge rearrange- 
ments originating in a variety of physical and/or chemical 
phenomena occurring between the electrode surface and the first 
organic layer can, in general, lead to appreciable ids^^'^^'^^'^^"^^. 
If available, for example, through a fully atomistic and quantum- 
mechanical treatment on a particular system, then they can, of 
course, be readily incorporated into the present model through 
equations (4) and (5). To study the impact of Wfgiec on Ae/h on 
more general grounds, however, we refer to an effective electrode 
work function 



<I>eff= Wfelec + ID, 



(6) 
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Figure 2 | Exemplary results of the iterative electrostatic model, (a) Calculated charge density p(z), and (b) electron potential energy —eV(z) within a 
d = 30nnn thick (^20 layer) pentacene filnn. Surface sensitive experimental methods (for example, UPS or Kelvin probe) only determine the total 
electron potential energy drop —eV(d) across the entire organic-semiconductor film, (c) —eV(z) for a range of different film thicknesses d up to 30 nm, 
highlighting that the series of experimentally determined -eV(d), shown as line connecting open circles, does not reproduce the band-bending 
— eV(z) within the organic semiconductor for any thickness. 




d (nm) 



> 

CD 



1 - 



10 



1 — 

20 



of (nm) 



30 




Figure 3 | Impact of effective substrate work function on band bending and barriers, (a) Total electron potential energy change —eV(d) across 
pentacene films of varying thickness d for a series of Fermi level energies ranging from 0.15 eV below E\_, the initial LUMO peak energy (top trace), to 
0.05 eV above £'h, the initial HOMO peak energy (bottom trace), (b) Data from a converted to injection barriers Ag for electrons (left axis) and Ah for holes 
(right axis). Dashed horizontal lines indicate respective values of zero and the vertical line highlights the cut visualized in c. (c) Injection barriers 
Ae (left axis) and Ah (right axis) as a function of effective substrate work function ^eff for a d = 10nm thick pentacene film with (•) and without (O) 
polaronic relaxation energies A + (for holes) and 1 ~ (for electrons) taken into account. The transitions between the regions with slope S = 0 and S = 1 mark 
the so-called pinning work functions, for LUMO pinning and <D+ for HOMO pinning. The lowest achievable injection barriers are not related to 
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in the following. Before charge equilibration, that is, with 
V{d) = 0, the initial barriers are then simply given by 
Ae = Oeff — EA and A^ = IE — Ogff, respectively. 

This now allows us to make the transition from relative 
potential energy changes for different values of (Fig. 3a) to 
absolute values of Ag and A^ for a series of effective electrode 
work functions, both of which are now quantities readily 
accessible through experiment. Figure 3b illustrates that changing 
Oeff indeed permits tuning the interfacial energy barriers/losses, 
but only within a certain range. For our model case of pentacene, 
the achievable minimum separation between and EA/IE 
rapidly converges to Ae = Ah = 0.4eV with increasing film 
thickness, where decreasing Ogff below ~3.0eV does not further 
reduce Ag and increasing Oeff beyond ~ 4.4 eV does not further 
reduce A^. A more convenient way of visualizing these data is to 
directly plot the energy barriers as a function of Ogff for a given 
film thickness, for example, d=lO nm as indicated by the vertical 
line in Fig. 3b. 

The obtained functional dependence (open symbols in Fig. 3c) 
is characterized by a slope of zero (S = 0) for low Ogff, a region 
with S = 1 for intermediate Ogff, and again S = 0 for high Ogff. The 
transitions between the different regions are marked by the so- 
called pinning work functions, O ~ = 3.0 eV for LUMO and O 
= 4.4eV for HOMO pinning with the latter matching the 
experimental value for a 10-nm film of pentacene^^. The term 
pinning refers to the observation that, for multilayer organic 
semiconductor films of device -relevant thickness, cannot cross 
the HOMO or LUMO level onsets (marked as horizontal lines) 
and, thus, energy barriers/losses cannot be completely eliminated; 
regardless of the choice for Ogff, the Fermi level appears to be 
'stuck' at specific energy separations ( ~ 0.4 eV for d—lO nm) well 
above/below the respective peak onsets. 

In the past, this observation has led to the postulation of 
material- specific, discrete, quantum-mechanical states at these 
energies, the so-called 'integer charge-transfer' (ICT) states, with 
ICT+ shifted up in energy from the HOMO onset and ICT~ 
shifted down in energy from the LUMO onset into the gap of the 
organic semiconductor^^. As already pointed out in refs 25,27,28, 
however, the fact that the pinning energies of E^ in the organic 
semiconductor gap emerge naturally (and quantitatively correct) 
from electrostatics alone: (i) lifts the necessity for invoking the 
existence of such states and (ii) invalidates the interpretation of 
their energy separation from the HOMO and LUMO onsets as 
the respective polaron-binding energies 1'^ and We 
explicitly illustrate this latter point in Fig. 3c for our Gaussian 
DOS {cf. ref. 28). For our considered model case, pentacene, the 
polaron-binding energy, that is, the energy stabilization of an 



excess charge in the organic semiconductor by geometric 
relaxation of the charged molecule itself (internal relaxation 
energy) and of the surrounding lattice (external relaxation 
energy), was found to be only =0.07eV (refs 57,58) with 

the latter being negligible compared with the former^^. This value 
is taken into account in equation (1) by shifting the HOMO and 
LUMO levels symmetrically towards each other, thus effectively 
reducing the gap of the organic semiconductor. The magnitude of 
the resulting changes (closed symbols in Fig. 3c) is, not 
unexpectedly, comparable to the absolute value of and, 
thus, its contribution to the observed minimum barriers is smaller 
than the typical experimental accuracy of ~0.1eV; polaronic 
effects will, therefore, no longer be considered in the following. 



Film thickness and Gaussian DOS. Briefly returning to Fig. 3b, it 
is evident that the relationship between Ag/h and Ogff changes 
with the total thickness of the organic film. In Fig. 4a, we show 
vertical cuts through the data presented in Fig. 3b not only for 
ii=10nm, as discussed in Fig. 3 c, but for an entire series of 
thicknesses. As a direct consequence of band bending within the 
organic film, the substrate work functions where LUMO/HOMO 
pinning sets in successively decrease/increase with decreasing film 
thickness down to 0~ = ~2.4eV and 0 + = ~4.7eV, respec- 
tively, for a monolayer of pentacene, which, again, compares 
favourably to experiment^^. Evidently, the thinner the organic 
semiconductor film, the lower Ag and A^ can be made by, 
respectively, decreasing or increasing Ogff until, in the extreme 
case of a monolayer, they vanish for Oeff=5.5eV and 
Oeff= 1.8 eV, respectively {cf. refs 30-33). 

The same could, arguably, be achieved also for thicker layers if 
the DOS tailing into the gap of the organic could be eliminated; it 
is the presence and the occupation of these gap states, after all, 
that prevents E^ from further approaching the HOMO/LUMO 
onsets. The origin of the gap states is commonly attributed to 
structural imperfection/defects^^'"*^'^^"^^. Their density and 
energy distribution can experimentally be determined only on a 
case-by-case basis^^"^^, as they depend on specifics of sample 
preparation and environment, that is, exact film structure and 
composition^^. To elucidate their impact on energy level 
alignment on more general grounds, we vary the width of the 
HOMO- and LUMO-related DOS while keeping it Gaussian- 
shaped, which has been shown to phenomenologically emulate 
also the impact of exponential tails . Leaving, at the same time, 
the energies of the peak maxima, E^ and £l> unchanged, trivially 
leads to the HOMO and LUMO onsets moving into the gap and, 
consequently, to getting pinned ever farther in the gap as well 
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Figure 4 | Impact of film thickness and width of Gaussian DOS. (a) Injection barriers Ag for electrons (left axis) and Ah for hole (right axis) as a 
function of effective substrate work function $eff fo'' pentacene filnns of thicknesses fronn d = 1.5 nm to 30 nm. (b) Injection barriers Ae (left axis) and Ah 
(right axis) as a function of ^eff for a d = 10nm thick pentacene film with the standard deviations of the Gaussian HOMO and LUMO DOS varied 
from (7h/l = 0.1 eV up to o-h/l = 0.4 eV while retaining a constant IE and EA. (c) Electron potential energy profiles — eV(z) at the highest effective substrate 
work function ^eff = 6.5eV for all o-h/l treated in b. 
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(not shown). A maybe less trivial result is presented in Fig. 4b. 
There, we keep the HOMO and LUMO onsets, which represent 
the physical material parameters IE and EA, respectively, at 
constant energy. Varying the width of the DOS then requires 
adjusting the energy separation between the peak maxima 
accordingly. Even under these conditions, however, reducing 
the Gaussian width of HOMO- and LUMO-related DOS leads to 
(i) the transition between the regimes of Fermi level pinning 
{S — 0) and the Schottky-Mott limit to become more 

abrupt, (ii) successively approaching the respective onsets in 
the pinning regimes and (iii) the spatial extent of the band- 
bending region to be reduced, as evidenced in Fig. 4c. The range 
of minimal energy barriers ( ~ 0.2 eV to ~ 0.6 eV) emerging from 
our model spans the entire range of values observed in 
experiment"*"*' , highlighting the decisive role of disorder in the 
organic semiconductor in determining interfacial energy level 
alignment^^. Note also that, for T>OK, finite barriers remain 
even in the limiting case cth/l^O of discrete molecular states, 
which is, in fact, the only case treated in refs 23,25. 

From band bending to interface dipole. Notably, the phenom- 
enology described so far, the characteristic shape in Figs 3c and 4 
in particular, is only ever observed for 'weakly' interacting 
electrode-organic interfaces. There, either through surface con- 
tamination or by design, the organic semiconductor is electro- 
nically decoupled from the metallic states by a thin, wide-gap 
buffer layer, which, in following refs 27,28, we did not explicitly 
take into account so far. As shown in the Supplementary 
Discussion, even doing so does not alter the fact that the inter- 
facial energy level alignment between organic semiconductor and 
metallic electrode always falls into one of three regimes, HOMO 
pinning, Schottky-Mott limit and LUMO pinning, characterized 
by slope parameters S = dAe/dOeff= — dA^/dOeff of (close to) 
zero, one and again zero, respectively (cf. Supplementary Figure 
1). If, however, the buffer layer is eliminated completely in 
experiments, which can be achieved by, for example, depositing 
the organic semiconductor directly onto atomically clean metal 
surfaces in ultrahigh vacuum, also fractional slope parameters 
between 0 and 1 are frequently observed^^'^^'^^'^^'^^. Despite 
having explored a good portion of the parameter space available 
in our model, we failed to generate non-integer values of S so far, 
whereas they are readily (and often quantitatively correct) 
predicted by the 'induced density of interface states' model^^'^^. 



This is surprising insofar, as this model is, in essence, quite 
similar to ours. Besides typically treating only molecular 
monolayers and assuming the organic DOS to be constant 
around for evaluating S, the main difference lies in the peak 
shape, which the authors take to be Lorentzian rather than 
Gaussian. This is entirely reasonable for the envisioned situation, 
where direct electronic coupling of the monolayer molecular 
levels to the metal leads to homogenous life- time broadening 
dominating over disorder- induced inhomogeneous (that is, 
Gaussian) broadening^^. Having already elucidated the impact 
of the latter in Fig. 4b, we now choose a Voigt peak function with 
successively increasing Lorentzian weight to represent the DOS of 
the first organic layer, while retaining pure Gaussians for all 
subsequent layers. The Oeff- dependence of the resulting injection 
barriers into 3. d= lOnm film of pentacene is shown in Fig. 5a, 
containing also the limiting, pure- Gaussian case already reported 
in Fig. 3c for comparison. 

Clearly, with increasing weight of the Lorentzian contribution, 
ever lower/higher substrate work functions are required to reach 
the same barriers for electrons/holes. More importantly, however, 
the Oeff-region of non-vanishing slope parameter expands, where 
S itself decreases from unity to a value of 0.42 for a Lorentz 
dominated (75%), down to 0.35 for an exclusively Lorentzian 
DOS, which brackets the experimental observed values of 0.37 
and 0.40 for pentacene^^'^^. 

This remarkable transition between phenomenologically dis- 
tinct energy level alignment regimes is accompanied by profound 
changes in the electron potential energy profile — eV{z) within 
the organic semiconductor film (Fig. 5b). On increasing the 
Lorentz weight in the first-layer DOS, the spatially extended 
band-bending region in the organic is successively compressed 
towards the immediate interface with the electrode until, in the 
case of a pure Lorentzian, it collapses into a potential energy step 
across only the very first molecular layer. The abruptness of this 
step, which is not a unique feature of the Lorentzian peak shape 
per se but simply a consequence of the (compared to the 
Gaussian) higher DOS around Ep, has led to it being assimilated 
into the interface dipole (as defined above) in the past. This can 
be justified by considering that the combined action of all 
interfacial effects reduces the work function that subsequently 
deposited molecules 'see' to, in our case, Oeff — ey((i) = 6.5 
eV — 2.1 eV = 4.4eV, which is, again, the HOMO-pinning work 
function O ^ of pentacene. Therefore, no further charge transfer 




<I>eff(eV) z(nm) <I>eff (eV) 

Figure 5 | Fractional slope parameters and band bending versus interface dipole. (a) Injection barriers Ae for electrons (left axis) and Ah for hole 
(right axis) as a function of effective substrate work function ^eff for a d = 10 nm pentacene filnn. The peak shape of the HOMO- and LUMO-derived DOS in 
the first layer is successively varied from pure Gaussian over Voigt to pure Lorentzian. (b) Electron potential energy —eV(z) within a d = 10nm pentacene 
film and a $eff of 5.3 eV. Changing the peak shape of the HOMO- and LUMO-derived DOS in the first layer from pure Gaussian over Voigt to pure 
Lorentzian leads to a collapse of the spatially extended band-bending region to an interface dipole across the first molecular layer, (c) Injection barriers 
Ag for electrons (left axis) and Ah for holes (right axis) as a function of effective substrate work function $eff for a d = 10 nm pentacene film. The DOS 
of the first layer is kept fixed to a Lorentzian peak shape of width 0.25 eV, whereas the standard deviations (Jh/l of the Gaussian peak shapes in all 
subsequent layers increase from 0.1 eV up to 0.4 eV. Although the range of linearity decreases succeedingly, the value of the slope parameter S<1 
remains unchanged. 
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between the monolayer-covered electrode and succeeding organic 
semiconductor layers occurs and, with having reached its 
minimum value of 0.4 eV already, no additional shifts in electron 
potential energy arise. Deviation from a strictly linear evolution of 
Ah/e is not observed until the Gaussian DOS of subsequent bulk 
layers again tails across for extreme values of Ogff, thereby 
causing residual band bending. The range of linearity is thus 
restricted by the amount of energetic disorder in the bulk 
(exemplified by crn/L-values ranging from 0.1 eV up to 0.4 eV in 
Fig. 5c), but the absolute value of the slope parameter S< 1 is not 
affected. This noteworthy result rationalizes the success of the 
induced density of interface state model in correctly reproducing 
the interfacial energy level alignment in the limiting case of 
atomically clean metals with intermediate work functions, despite 
it considering only monolayer coverage. 

Discussion 

To summarize, we present the results of a single, versatile 
electrostatic model, from which all observed regimes of energy 
level alignment at metal/ (insulator/) organic semiconductor 
interfaces emerge naturally. Although designed to accommodate 
unprecedented detail and complexity, it contains prior work as 
limiting cases, notably, that of Ley et al?^, which is limited to 
producing either VL alignment (S=l) or Fermi level pinning 
{S — 0) as a consequence of considering only discrete molecular 
states; that of Blakesley et a/.^^'^^, which does not explicitly 
include insulating buffer layers and, thus, overestimates 
interfacial charge densities, as shown in the Supplementary 
Discussion; and that of Vazquez et al?^, which is limited to a 
single molecular layer with Lorentzian DOS and, thus, fails to 
predict band bending. In providing a unified picture of the 
underlying physical mechanisms, our model confirms that the 
DOS of the organic semiconductor is the central, decisive 
quantity, responsible for all experimentally encountered 
phenomena and that controlling it is, therefore, key to 
minimizing charge carrier injection barriers and/or extraction 
losses beyond what is possible by adjusting the electrode work 
function alone. 

Using high (for holes) and low (for electrons) work-function 
electrodes is, of course, still a necessary prerequisite. Our results 
illustrate, however, that there are intrinsic limits to this strategy 
and that, to overcome these limits, the amount of electronic states 
in the energy gap of the organic semiconductor (as well as their 
energetic separation from its band edges) must be reduced. 
Besides eliminating defects within the organic semiconductor, our 
model suggests electronic decoupling from the metallic electrode 
to be a critical component in achieving this, thereby providing a 
theoretically well-founded justification for the wide-spread use of 
thin passivating buffer layers. As a further beneficial consequence 
of 'cleaning out the gap', spatially extended band-bending regions 
are predicted to develop within the organic semiconductor, 
which, in contrast to sharply localized interface dipoles, further 
lower effective injection barriers in OLEDs under applied bias 
voltage and facilitate the efficient collection of photo -generated 
charge carriers in OPVCs by imposing a drift direction. Our 
results thus provide a solid theoretical basis for past and present 
efforts in optimizing the interfacial energy level alignment at 
organic semiconductor/electrode contacts and will, hopefully, 
inspire future work. 
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